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ABSTRACT 

Compared with the cold dark matter (CDM) model, in the warm dark mat- 
ter (WDM) model formation of small scale structure is suppressed. It is often 
thought that this would delay the reionization of the intergalactic medium (IGM), 
as the star formation rate during the epoch of reionization (EOR) would be low- 
ered. However, during the later stage of the EOR, a large portion of the ionizing 
photons are consumed by recombination inside the minihalos, where the gas has 
higher density and recombination rates than the gas in the IGM. The suppression 
of small scale structure would therefore reduce the recombination rate, and could 
potentially shorten the reionization process. This effect is investigated here by 
using the analytical "bubble model" of reionization. We find that in some cases, 
though the initiation of the EOR is delayed in the WDM model, its completion 
could be even earlier than the CDM case, but the effect is generally small. We 
obtain limits on the WDM particles mass for different reionization redshifts. 

Subject headings: dark matter — large-scale structure of universe 



Introduction 



The nature of the dark matter is presently unknown. The cold dark matter (CDM) 
model has been very successful in explaining the observed properties of large scale structures 
of the Universe, but on small scales there are still some discrepancies, as the abundance 



of satellite galaxies falls far short of th e number of subhalos predicted (Klypin et al. 



Moore et al. 



19991 ; ISpringel et al.ll2008l ; but see also ISimon fc Gehal 120071 and iPrimack 



999; 



2009 



for a different view). An interesting alternative is the warm dark matter (WDM) scenario 



where the dark matter particle has a smaller (a few keV) mass (IBlumenthal et al.l 11982 
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). As the growth 



of structure below the free-streaming scale of the WDM was quenched, such a model predicts 
much fewer small halos, but on large scales its prediction is similar to the CDM cas e and 



could still match the l arge scale structure observations very well (IBode et al. 
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The abundance of small halos may also greatly affect the reionization of the Universe. 
As the stars and galaxies form at the end of the cosmic dark age, the star light ionizes the 
intergalactic medium (IGM), causing the later to ionize. Eventual ly, all of hydrogen gas 



in the IGM is i onized, this is the so called (hydrogen) reionization (IBarkana fc Loebl 12001 
Fan et al.l 120061 ) . In the hierarchical structure formation scenario, most of these ionizing 



photons would come from dwarf galaxies in small halos. In recent years, the Gunn-Peterson 
trough in the Lyman alpha absorption spectrum of high redshift quasars has been observed, 



Becker et al. 


200ll: 


Fan et al. 


2002) 



anisotropy observation is another important source of information about EOR (IZaldarriaga 



19971 : iKogut et al.ll2003l : ISpergel et al.l 120071 ). The recent WMAP anal ysis gives the redshift 



of completion of reionization as z re = 10.5 ± 1.2 (ILarson et al.l 1201 ll ) . This result can be 
used to constrain the properties of dark matters. In the WDM model, halos of a given mass 
would typically form later than in the CDM model, so the initiation of the EOR would be 
delayed. Also, as the formation of small halos which are most common during the EOR are 
suppressed, fewer ionizing photons are produced, and it was believed t hat this would delay 



the completion of the reionization significantly ( IBarkana et al.l 1200 ll ; lYoshida et al.l 12003 
Somerville et al.l 120031 ). 



It should be noted that in some WDM models, the WDM particle may be able to de 



cay, which ionizes and heats up the IGM (lAbazajian et al.ll200ll ; iMapelli fc Ferraral 12005 



Mapelli et al.ll2006l ). This (and the annihilation of some CDM particles) in itself is an inter- 



; Chuzhov 


2008; 


Yuan et al. 


2010: 



20101). The effect on reionization is 



fairly complicated, for example, the free electrons produced in the dark age could help cat- 
alyze the formation of molecul e hydrogen, making it possible to form first stars and start 



reion ization early in this case (IBiermann fc Kusenko 



2006; Kusenko 2007; Stasielak et al. 



2007h. though the magnitude and significance of this effect is highly uncertain( ]Ripamonti et al.l 



2007al lbh. In the present paper, we shall limit ourselves to the passive WDM model, and 
shall not consider such effects. 
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However, the time it took the reionization process to complete depends not only on the 
supply rate but also on the consumption rate of the ionizing photons. After being ionized, the 
atoms of the IGM gas will recombine, and to keep the gas ionized more ionizing photons have 
to be consumed. As the recombination rate is proportional to the density squared, at the 
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2010). In the 



WDM scenario, the number of minihalos are much fewer, so the global recombination rate 
could be significantly lower. When this is considered, it is not obvious whether the WDM 
would delay or advance the completion of the EOR. 



We use the bubble model of reionization fjFurlanetto et al.l l2004al ) to investigate the 
reionization process in the WDM scenario, take into account both the reduction in photon 
production rate and the consumption rate due to the suppression on halo formation in WDM 
models. Here we adopt the WMAP 5 years cosmology parameters (the WMAP 7 years pa- 
rameters are almost id entical): (ft m , Q A , ft 6 , h, a 8 , n s ) = (0.274,0.726,0.0456,0.705,0.812,0.95) 
( iKomatsu et al.l 120091 ) . Note that for the WDM model, this set of parameters may not be 
the best fit values, but to illustrate the physical effect of different dark matter mass, we have 
used the same parameters for all models. 



2. Method of Calculation 



We calcul ate the halo mass functio n in the WDM model analytically by following the 
prescription of ISmith fc Markovid ( 120111 ). In the WDM model, the free-streaming comoving 
scale is given by 
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and the corresponding mass is Mf s = 47r/3(A/< j /2) 3 p m . The halo mass function is then given 
by 

dn , , , , 1 



dM 



{M,z) 



1 + erf 
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PS 



where aing iu = 0-5, and [jjjj] rc can be calculated wi th the usual Press-Schechter (PS) pre- 



scription ( IPress fc Schechterl Il974l; iBond et al 



2001 



1991 ). with the m atter power spectrum for 



Viel et all 120051 ) 



the WDM case given by the fit (IBode et al 

PwMk) = Pcdm(A;){[1 + {akf^^} 2 



(3) 
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where /i = 1.12, and 



a 



0.049 



V keV J 



-i.il 



^WDM 

0.25 



0.15 



/l 

0/7 



1.22 



h~ l Mpc. 



(4) 



The halo mass functions for the CDM model and the WDM model with different particle 
masses at redshift 15 are plotted in Fig. [TJ We see the abundance of dark matter halos on 
small scales are significantly suppressed, while on large scales the halo abundance are almost 
the same as the CDM model. On the same plot, we also mark the mass range of what we call 
the "minihalos", i.e. halos which are sufficiently massive t o accrete the gas, b ut not massive 
enough for the gas to cool by atomic H and form galaxies tlCiardi et al.ll2006l ). 
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Fig. 1. — The halo masses function for CDM (solid line) and WDM with m\yDM = 10, 2 keV 
(dashed, dashed-dotted) respectively at z=15. We marked the minihalos mass range by two 
thin dashed vertical lines. 



We now consider the reionization process in the WDM model. At the end of the dark age, 
dark matter halos form from primordial overdense perturbations, and in halos of sufficient 
masses, the gas cooled radiatively and contracted to dense clumps which eventually form 
the first stars. At present, it is still unclear whether the first stars formed in halos which 
were cooled by the trace amount of molecule hydrogen in the gas (T v ; r >~ 10 3 K), or only in 
more massive halos which were cooled by atomic hydrogen (T vir > 10 4 K) f lBarkana fc Loeb 
200ll ; iBromm fc Yoshidal l201ll ). In either case, the stars and galaxies tend to form more 



abundantly in overdense regions, so the ionized regions first appeared as "bubbles" in such 
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regions. As more and more bubbles appeared and grew in size, eventually the ionized regions 
overlapped, thus completing the reionization. 

Inspired by numerical simulation res ults, a popular analytica l model of the process cal led 
the "bubble model" has been developed ( iZaldarriaga et al.l 12004 ; iFurlanetto et al.ll2004aU bl). 
According to this model, the reionization proceeds by first forming ionized bubbles, the 
number and size of these bubbles continue to grow till they overlap with each other, thus 
complete the reionization process. The formation of the ionized bubbles is determined by 
the condition that the number of ionizing photons produced within the given region exceeds 
the total number of atoms to be ionized in the same region. If we assume that each collapsed 
baryon on average contribute ( ionizing photons, with each photon ionizes one atom once, and 
the average number of recombination per atom during that time is n rec , then the condition 
for the region to be ionized at the given redshift can be written as 



Cfcoll > 1 



n r 



(5) 



where f co u is the fraction of baryons collapsed into sta r -form ing: halos. This treatment of the 
effect of recombination differs from IFurlanetto fc Ohl (120051 ). They considered the limiting 
case, where the ionizing photon production is required to counteract the recombination at 
that instant. Our model in Eq. dSJ) is more in line with the original bubble model, in which 
the total integrated number of ionizing photons is considered. 

In such a model, one assumes that halos with the virial masses above a certain threshold 
value could form stars. The collapse fra c tion can then be calculated with the extended 
Press-Schechter method faond et al.lll99ll : lLacev fc Cold 1199a IMo fc White! Il996h . During 
the early EOR, the Pop III stars formed in the molecule hydrogen cooled halos may have 
played an important role, these halos are much less massive than the typical star-forming 
halos during the EOR, which are cooled by atomic hydrogen or at somewhat later time, 
by metals. To account for this, we may split the contribution into the part from molecule 
hydrogen cooled halos and the part from more massive halos. We denote the first type by 
the subscript "mol", and without danger of confusion, the second type is denoted without 
any subscript, then the condition Eq.flSJ) can be rewritten as 



(mol f mol + Cf > 1 



(6) 



A molecule hydrogen cooled halo should have a virial temperature of at least 10 3 K for 
this mechanism to work. Furthermore, the formation of the molecules is strongly modulated 
by the Lyman- Werner (LW) rad iation background . For the molecules to form, the mass of 
the halo should be greater than (ITrenti fc Stiavellil 120091 ) 

.1 + *, 
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The mass threshold is then given by 

M mol (z) = max[M^ r (10 3 K, z), M H2>cool (J 2 i, z)\. 
We model the LW intensity as J = J%\ x 10 _21 ergs _1 cm _2 Hz _1 sr _1 with 



Mz) = 10- 3 \^f) + 0.28 (^j (1 + zff coll {z), (9) 

The first term represents a slowly-evolving component due to Pop III stars whose formation 
is self-regulated, while the second term comes from the more massive, atomically cooled 
halos. Strictly speaking, for the WDM model the formation of Pop III stars is suppressed, 
so the coefficient of the first term would be smaller, but this change would not affect the 
final result, as the contribution of the few stars in WDM model is very small anyway The 
molecule cooled halo collapse fraction is then calculated with 

i r rM vlr (W* K,z') j2 

Uoi = — dz' M———dM. (10) 

Pm J JM mol {z<) dMdz' 

The value of ( and ( mo i depends on a number of factors, e.g. the fraction of baryons 
in the halo which ended up in stars, the energy released by the star during its lifetime, 
and the fraction of the photons which escaped to the IGM from the cloud surrounding the 
star. There are curren tly large uncertainties on this parameter. Here we adopt ( = 40 



(IBarkana fc Loebll200ll ) and ( mo i = 50, the l atter parameter value is in agreement with the 



properties of massive metal-free stars given in lSchaererl (120021 ). With these parameter values, 
J21 ~ 10 -3 during the Pop III stars dominated stage in the CDM model, in w ell agreement 



with the value widely adopted in analytical and numerical stu dies of first stars f|Wise fc Abel 
2007 ). while at z=10 J21 ~ 7, very close to the value given in Trenti &: Stiavelli f 20091 ). 



At high redshifts, the clustering in the IGM is relatively low, we assume that only one 
photon is needed to ionize an IGM atom, while for the atoms in the minihalo the number of 
recombinations is given by 

1 fMup(z) dn 

n rec , MH = — / iM—dM. (11) 

Here £ is the average number of recombinations per atom of the minihalo, Mj is the Jeans 
mass. The upper limit of the integration M up is the minimum mass of halos that could host 
radiation sources, i.e., Pop III stars or galaxies, that is mm[M mo i(z), M vir (10 4 K,z)]. The 
recombination inside the galaxies should be relegated to the net photon production number 
( and should not be included again here. Earlier analytical estimates typically gave £ ~ 10 2 
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( jHaiman et al.ll200ll ). in which case the recombination of the minihalos would consume the 
majority of photons at the late stage of EOR. However, numerical simulation of the ionizing 
front passing through minihalos shows that the number of actual recombinations may be far 
smaller, e.g., at z=15, £ < 8 for halos with mass below 10 7 M Q and irradiated by typical 



ionizing flux (IShapiro et al.l 12004 ; llliev et al.l 120051 ). As the halo is ionized, it is heated and 
the gas expands, causing the density to decrease and recombination rate quickly lowered, 
and thus the recomb ination rate is much lowe r than the earlier estimates. We adopt the £ 
value derived from Shapiro et al. ( 2004 ) and Iliev et al. ( 2005 ). 



In Eq. 06]), for a bubble with mass m, f and f mo i are both functions of the linear 
overdensity of this bubble. By solving the equation Cmaifmoi(z, S x ) + Cf( z J S x ) = 1 + n rec , a 
critical overdensity S x (m,z) is determined. This is the barrier in the excursion set, re gions 
with linear overdensity above this barrier should be ionized ( IFurlanetto et al.ll2004al ). We 
obtained this barrier by solving the above equation with numerical iterations, and found that 
it is still well approximated by a linear functi on of the squared varian ce of density fluctuation 



a 2 (m) (the following process is the same as in IFurlanetto et al.ll2004al ): S : 



m, z) 



B(m, z) 



B + Bicr 2 (m), where B = 5 x (m — > oo, z) and B 1 = |%(m — > oo, z). With this linear barrier 
the bubble mass function is expressed analytically by 



dn b 
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2 p r , 



7r m z 
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exp 
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(12) 



where rib is the number density of bubbles. Finally, the volume filling factor of all bubbles 
Qv is calculated directly by integrating over the bubble mass function: 



V(m)—^dm. 
dm 



(13) 



3. Results and Discussions 

We plot the redshift evolution of the ionization volume filling factor in Fig. [2j Here we 
considered four different cases: the reionization in the ACDM model with and without mini- 
halo recombinations, and the reionization in AWDM model with minihalo recombinations 
for mwDM = 10 keV and 2 keV. 

First we note that as illustrated in the ACDM case, in the absence of minihalos, the 
reionization would be completed earlier by as much as Az = 1, this shows how much impact 



1 Our definitio n of £ differs by one from 
Iliev et al.l |2005l) . 



Iliev et al. 



(|2005|) : £ = £iii 0V — 1, where £ni ov is the value given in 
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the minihalo rec o mbina tions could have on the reionization. This is obtained with £ derived 
from llliev et al.l (120051) . If we adopt the higher values of £ given in the earlier literature 
(IHaiman et al.ll200l[ ). the impact would be even stronger. 




Fig. 2. — The redshift evolution of the ionized bubble volume filling factor in different models. 

In the AWDM case, the reionization starts later than in the ACDM cases, as the first 
stars would form later with the small scale power suppressed in these models. This difference 
is most obvious at z > 20, when the first stars just begin to form in large numbers. As 
small scale powers are suppressed in the AWDM model, there are very few halos in which 



the g as could be cooled by molecule hydrogen and form Pop III stars (10 'Shea fc Norman 



20061 ). Instead, in AWDM most first stars only form in atomically cooled halos, this caused 



significant difference in the initiation of reionization. 

However, it is generally believed that the reionization is due primarily not to the stars 
formed in the H 2 cooled halos, but to stars formed in the more massive halos at lower 
redshifts. For the more massive halos, the difference between the lOkeV WDM and the 
CDM is not so large. In fact, from the figure we can see that in this case, the bubble filling 
factor of the lOkeV WDM model catches up with the CDM model at z ~ 13, and later it 
even exceeds that of the CDM model, thus the reionization is actually completed earlier than 
in the CDM case. This is because at this stage, the ionizing photons are produced mainly 
in the more massive halos, which is about equally abundant in the 10 keV WDM and CDM 
cases, while the suppression of minihalos reduced the global recombination rate in the WDM 
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model, making the bubbles overlap earlier. Thus we observe the interesting result that the 
Universe is reionized earlier in the WDM model. This effect is relatively small compared 
with the current theoretical and observational uncertainties, especially the mean number of 
ionizing photons produced by a halo of given mass. However, it could potentially be useful 
if these uncertainties are greatly reduced by improvements on observations and theoretical 
modelings of reionization process. 

For models with still lower WDM masses, the more massive halos are begin to be 
affected. In the 2 keV WDM model, the star formation associated with the more massive 
halos are also suppressed, and the reionization is delayed in that model. The reduction of 
photon supply and consumption would be balanced somewhere between these two cases, the 
exact value would however be dependent on the parameters we adopt, especially the value of 
the parameters ( and £. The values we adopted here are plausible and non-extreme, with the 
reionization happen at z ~ 10, in agreement with the WMAP constraint. Nevertheless, there 
are large uncertainties in these parameters. In particular, if the £ value for the minihalos is 
greater, as was used in some earlier papers, the effect would be still stronger. 

We plot the redshift of completion of reionization as a function of WDM particle mass 
in Fig. |3l For reference, the corresponding cutoff mass scale M cut , where the mass function 
dn/dM is suppressed by a factor of e when compared with the ACDM model, is also plotted 
on the upper abscissa. We see that below 11 keV, the reionization redshift is very sensitive to 
the WDM particle mass. With jtiwdm increases, the reionization redshift also raises quickly, 
as was usually assumed for WDM particles. However, a peak of reionization redshift of 
z re 9 is reached at m^DM ~ llkeV. Above this mass, the reionization redshift begin 
to decrease slowly, as the suppression on number of star-forming halos becomes relatively 
insignificant, while the suppression of minihalos reduced the global recombination rate. As a 
comparison, we also plot z re in the CDM model with and without minihalo recombinations 
by the filled and open circles respectively. Quasar absorp tion line studies show that the 
Universe had been reionized at least ~ 6 (IFan et al.ll2006l ). This gives the constraint that 
mwDM > 1.3 keV. 

This limit on WDM mass may b e compared with r ecent WDM mass limits obtained 
from other observations. For example, iNaravanan et al.l (120001 ) obtained m^DM > 0.75 kev 
from the Lya forest observations, while 



VieletaL 



(200 5j) obta i ned a lower limit of 0.55 kev 



from CMB (WMAP) and the Lya forest data. iBarkana et all (120011) gave m W nM > 1 .2 keV 
with the requirement z re > 5.8 for their fiducial model, while iPolisensky &: Ricottil (120111 ) 
obtained a constraint of m^DM > 2.3 keV from the number of Milky Way satellites. 

In the above, w e have adopt the same value of ( as used in the fiducial model of 
Barkana et al.l (l200l[ ). which was derived from the observations of z ~ 3 — 4 and present- 
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Fig. 3. — The redshift of the completion of reionization as a function of wwdm (the solid 
line). On the right side of the figure, the open circle represents the result for the CDM model 
without minihalos, while the filled circle represents the result for the CDM model when the 
effect of minihalo recombination is taken into account. For comparison, the mass scale M cut , 
where the mass function dn/dM is 1/e of the LCDM value, is also marked on the upper 
abscissa of the figure. The shadowed region corresponds to z re < 6.0, already excluded by 
the observations of Gunn-Peterson trough in z ~ 6 quasar absorption line spectra. The gray 
region above is the constraint from the WMAP, i.e. z Te = 10.5 ± 1.2. Here we adopted the 
fiducial value ( = 40. 



day galaxies. However, there are still large uncertainties on the properties of sources in 
the epoch of reionization. If these sources are stars with metallicit y Z = 5 x 10 -4 Z@, 
and if their initial mass function is the Salpeter form ( jSalpeterl Il955l ) with the mass range 
1 < M < 100M o , for t he starburst m odel the number of ionizing photons produced per 
stellar atom is ~ 13000 (jSchaererll2003l ). Assuming a star formation efficiency of 0.05 and 
escape fraction 0.5, ( could be as high as ~ 300. Considering these uncertainties, we also 
calculated the reionization redshift for different ( values. In Fig. HI we plot the contours of 
z re with both ?7i wdm and ( varying. We see WDM models with particles mass below 0.5 keV 
have already been excluded, otherwise the reionization could not be completed before redshift 
6 for reasonable values of (. On the other hand, for particles with m > 8keV, ( should be 
less than ^180, otherwise the completion of reionization would be too early, in conflict with 
the WMAP observations (ILarson et al.ll201ll ). 
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Fig. 4. — The contour map of the redshift when reionization is completed in the m\yDM-C 
plane, we label three lines, z re = 6, which is the lower limit of reionization redshift given by 
Gunn-Peterson trough observations, and z re = 10.5 ± 1.2, which are from WMAP. 



Besides the halo abundance, the distribution of gas within the halo is also affected 
by the replacement of CDM with WDM. The collapse of halos in the WDM model is 



generally later than in the CDM mode 



hen ce the halo concentration is typically smaller 



(IBarkana et al.ll200ll ; ISmith fc Markovid 120111). as the dark matter halos formed later would 



reach s maller average densi t ies (jNavarro et al.lll997l ). In addition to the delayed collapse 
redshift, ISmith &: Markovid (120111 ) pointed out that due to the relic velocities of WDM par- 
ticles, the core of the halo would be smoother in the WDM case. These effects could in 
principle reduce the photon consumption rate of the minihalos. We have checked this effect 
by calculation and found that it should only have a very slight influence on our results. 
For example, in the ttiwdm = 10 keV case, for a halo with 10 7 M Q at z — 10, the density 
profile of the halo within 4 x 10~ 2 R vir would be flattened, but outside this core radius the 
density profile is hardly changed. The gas density profile changed even less, since even in 
the CDM model the baryonic gas has pressure and therefore a more smooth distribution. 
The larger minihalos contribute more to photon consumption in the WDM model, because 
the abundance of smaller ones is reduced significantly. In the end, we find that the change 
in density profile does not significantly affect the total recombination rate of minihalos. 



In the present calculation, we have used the bubble model to treat the reionization 
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process. The bubble model is of course only an approximate model, though it does repro - 



duce more elaborated simulations (jZahn et al.ll2007l ; iMcQuinn et al.ll2007l ; IZahn et al.ll201ll ). 
Furthermore, the bubble model only provides a distribution of bubbles at a given redshift, 
it does not tell us how each individual bubble have grown, so we have made only an static 
treatment of the recombination process: we calculated the volume of ionized bubble at a 
certain redshift, then counted the number of minihalos in this region at that time, and 
calculated the number of recombinations for these minihalos until they are all photoevap- 
orated. This treatment might slightly overestimate the number of recombinations, because 
some of the bubbles are formed by the growth and merger of smaller bubbles which were 
formed at earlier time, and in these regions the minihalo formation had been suppressed. To 
check whether the neglecting of bubble growth history would alter our results, we make a 
simplifying ansatz of bubble growth. We assume that for a given bubble at redshift z, the 
different parts of its volume were acquired (i.e. first ionized) at different redshifts z' > z, and 



the contribution at z' is proportional to 



dz 



With this assumption on bubble growth 



history, the suppression on minihalo formation in previously ionized region can be taken into 
account, and the average number of recombinations within the bubbles is given by 



1 



n rec,MH 



fcoll\z) 



n rec,MH(z')- 



dfcoii(z'] 



dz' 



dz'. 



where n reCt MH on the R.H.S is given by Eq. f fTTT) . We then calculate the evolution of Qy with 
this new recombination number. However, we find that the resulting difference is very small. 
This is not surprising, since increases rapidly as redshift decreases, in the calculation of 



n r ec,MH{Z), n rec MH{ z 



contributes most to the integration. 



Finally, iGao fc Theunsl (120071 ) suggested that in the WDM models, unlike the case of 
the CDM models, the first stars do not necessarily form in halos, but may instead form in 
filaments. They also argued that with this new way of star formation, the first stars may have 
smaller typical mass (~ M ), and the global star formation rate could be even higher. This 
also raises the interesting possibility that the reionization could occur earlier in the WDM 
model than in the CDM model. However, whether this new formation mechanism would 
indeed work as they proposed is still not completely clear, as it is very difficult to model the 
star formation process in sufficient resolution at present to really check the outcome of the 
filament star formation process. Here, we have maintained the standard view, and assumed 
that the star s formed only in DM h alos. We note that if the stars indeed form in the way 
suggested by IGao Sz Theunsl (120071 ). the reionization in the WDM model would occur at 
even higher redshifts, and there is not anything incompatible with the recombination effects 
we discussed. 



In summary, we find that in WDM models the completion of reionization is not always 
delayed, in some cases, when the WDM mass is not too low, it could even be advanced 
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due to the reduction of recombination rates. This effect is relatively small compared to the 
uncertainties in observation and theoretical models at present. We also find that for ( = 40, 
to be consistent with the observations of Gunn-Peterson trough in z ~ 6 quasar spectra, 
the mass of dark matter particles should be higher than 1.3 keV. However, if more ionizing 
photons were produced and escaped into the IGM, the dark matter could be warmer, i.e 
have smaller mass. However, for dark matter particles that are less than 0.5 keV, ( > 500 
is needed. This would only be possible if the reionization photons are contributed primarily 
by metal-free or massive extremely metal-poor stars. 
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